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Abstract In this paper, the rms-flux (root mean 
square-flux) relation along the Z-track of the bright 
Z-Source Cyg X-2 is analyzed using the observational 
data of Rossi X-ray Timing Explorer (RXTE). Three 
types of rms-flux relations, i.e. positive, negative, and 
'arch'-like correlations are found in different branches. 
The rms is positively correlated with flux in normal 
branch (NB), but anti-correlated in the vertical hori- 
zontal branch (VHB). The rms-flux relation shows an 
'arch'-like shape in the horizontal branch (HB). We also 
try to explain this phenomenon using existing models. 

Keywords stars: individual (Cyg X-2)-stars: neutron-J 
X-rays: stars 



1 Introduction 

Aperiodic X-ray variability is a common property 
in X-ray bi naries (XRBs) an d active galactic nu- 
clei (AGN) (Ivan der Klisl Il994 IVaughan et ail 12003 : 



McHardv et al.ll2004T ). Many works have been done to 
investigate the aperiodic X-ray variability. For exam- 
ple, a linear correlation between rms (see Equation 1) 
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an d flux of the X-ray in AGN and XRBs is reported 
bv lUttlev and McHardvl (|200l[h This relation is found 
in both neutron st ar XRBs (NSXRBs) and b l ack hole 



201(j 



XRB s fBHXRBs) (lUttlev and M cHardv 120011 : iLi et al 



Gleissner et ah 2004; 



Heil et al. 2011). Accord- 



ingly, the linear relation is supposed to be a fundamen- 
tal pro perty of accretion flows in compact objects. Re- 
cently, (Heil and Vaughanl (|2010l ) report that this linear 
relation also exists in a Ultraluminous X-ray source, 
NGC 5408 X-l, which further strengthens the above 
supposition. 

The rms-flux relation is very useful in evaluat- 
ing X-ray variability models. Considering that the 
shot-noise model predicts a stationary power spec- 
trum and can't pro d uce a linear rms-flux relation, 
Uttlev and McHardvl (|200lh pointed out that this 



model is not correct. Furthermore, lUttlev! ([20041 ) 
showed that in the accreting millisecond pulsar SAX 
J1808.4-3658 is coupled with the 401-Hz pulsation. 
This indicates that the linear relation is not produced in 
the corona but in the accretion flow onto the magnetic 
caps of the neutron star (NS), which means that the 
linear relation does not favor the coronal flare model 
for the X-ray variability neither. Based on the in- 
vestigation of the rms-flux line ar relation on different 
time-scales. lUttlev et al.l (|2005f l ruled out several X-ray 
variability models except the perturbed ac creting-flow 
model (e.g. density waves, iMisral rt2000h of Lvubarskii 
( 1997t l . On the other hand, IZhaiiJ |2007|) found a 
similar rms-flux relation in the Sun and 7-ray bursts, 
suggesting that the accreting-flow model may be also 
wrong for the X-ray emission since the Sun does not 
have accretion flows. In conclusion, the rms-flux rela- 
tion is very important in distinguishing X-ray variabil- 
it y models. 

Gleissner et al. ( 2004( l investigated the rms-flux rela- 



tion in all known spectral states of a BHXRB, Cyg X-l, 
and found that the linear relation exists in all spectral 
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states although their po wer density spectra (P DS) show 
rather different shapes ( Gleissner et al. 20041 ). This re- 
sult implies that the rms-flux relation is a more univer- 
sal characteristic than PDS shape. 

Cyg X -2 is a bright, persistent l ow mass XRB 
CLM XB1 dBradt and McClintockl [l983h . It contains a 
NS ( Kahn and Grindlavl 1984 ). with an optically mea- 
sured massofj^l^JSiO^M^and distance of 7.2±1.1 
kpc ( Orosz and Kuulkers 19991) . The orbital perio d 
is ~9.8 days (|Cowlev et al.lll979l; ICasares et ai1ll998l) . 
For comparing the differences between a BHXRB and 
a NSXRB, we will investigate the rms-flux relations of 
different states in Cyg X-2 to check out whether the 
linear relation is valid or not throughout all the states 
of a NSXRB. 

Cyg X-2 is a Z source according to its track through 
differen t spectral states on the color-color diagram 
(CCD) dHasinger and van der Klisll989llHasingerll990l; | 



van der Klisl 12000). In general, a Z source has three 
branches, i.e. the horizontal branch (HB), normal 
branch (NB) and flaring branch (FB). The HB usu- 
ally locates at the top of the CCD, followed by the NB 
crossing the diagonal, and the so-called FB lies at the 
bottom of the CCD. Usually, a Z source can trace out a 
Z-track in few hours, but the Z-track will shift in weeks 
or months. Sometimes, the HB may have two portions: 
one is horizontal and the other is vertical. We call 
this vertical portion vertical horizontal branch (VHB ) 
( Shirevl[l99lilKuulkers et al.lll999t Ishirev et al.lll999l ). 
During the evolution of a Z-source, it moves continu- 
ously on the Z-track, du e to different mass accretion 
rate (|van der Klisl |2000|) . From simultaneous multi- 
band observations, it is thought that the mass accretion 
rate reaches a minimum at the left end of the HB, in- 
creases as the source travels from HB, throughout NB, 
to F B, and reaches its maximum at the right end of th e 
FB dHasinger and van der Klisl Il989i: lHasingerl Il990l ). 
However, Church et al. ( 20061) ha s a different opinion 
about the mass accretion rate (see Church et al. ( 2006f) 
and reference therein). 

Cyg X-2 varies on time-sca l es from milliseconds t o 
months ( Kuulkers et al.l 1996 ; Wiinands et al. 199G). 
Several kinds of quasi-periodic oscillations (QPOs) are 
detected in Cyg X-2 and its rapid variability is also 
found to be related with its position in the CCD. The 
QPOs that appear in the HB and the upper NB are 
called HB oscillations with the centroid frequency vary- 
ing from ~ 10 Hz to ~ 70 Hz. The NB oscillations 
locate in the lower NB with the centroid frequency at 
5-7 Hz (Hasinger and van der Klislll989t van der Klis 



2000; IPiraino et al.ll2002l: Ivan der Klisl |2004) . In addi 



tion, some kHz QPOs are observed in Cyg X-2 when the 
sour ce moves from the left end of the HB to the upper 
NB (|Wiinands et al.l (|1998[ ) and references therein). 



Section 2 is data reduction. The rms-flux relations 
are given in Section 3. Discussions are presented in 
Section 4. 



2 Observations and data analysis 

The observational data we used here are from RXTE 
archive, and are reduced by standard Heasoft package 
released by HE AS ARC. 

2.1 Observations 

In order to investigate the rms-flux relation in differ- 
ent branches, we require the exposure time as long as 
possible. Taking into account that the time interval 
Cyg X-2 stays in only one branch is ~ 10 000 s, we use 
only those observations whose duration is longer than 
10 000 s. Moreover, we discard those observations with 
different number of PCUs switched on during the ob- 
servation (see The RXTE Cook Book) and also those 
observations locate at the FB or at the lower part of 
NB, since their PDS always have negative Power Den- 
sity in the frequency range we are interested in. We list 
the observation IDs of the satisfactory observations in 
Table [B 

The first 7 observations belong to Epoch II and the 
other are during Epoch III. Their durations are longer 
than lOks. In ObsIDs 10065-02-01-000 and 10065-02- 
01-001, the event mode data, E_125us_64M_0_ls, are 
used for investigation while the binned mode data, 
B_2ms_16A_0_35_Q, are used in the other 14 observa- 
tions. We also list the location by their CCD position 
and PDS. 

The CCDs of each observation are made by using 
the Standard 2 mode data with time resolution 16 s. 
We define the soft color as the count-rate ratio between 
4.7-7.5 kcV and 2.0-4.7 keV, and the hard color be- 
tween 11 2-18.8 keV and 7.5-11.2 kcV during Epoch II 
(see also Lei et al. ( 20081 ) V Because the PCU gains are 
different for different Epoches (see The RXTE Cook 
Book) , the soft and hard colors of Epoch III are defined 
a little different from that of Epoch II, i.e., the count- 
rate ratio between 4.1-6.2 keV and 2.0-4.1 keV, and 
between 9.4-15.9 keV and 6.2-9.4 keV for soft and hard 
colors respectively. However, the difference is corrected 
using the Crab observations. All the above count rates 
are background-subtracted. The PDS of all the obser- 
vations are extracted from the Binned Mode or Event 
Encoded mode (see Table [TJ to locate their positions in 
CCD. The branches of those observations are listed in 
Tabic [TJ 
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2.2 rms calculation 



3 results 



The Binned Mode or Event Encoded mode data are 
also used to obtain the rms-flux relations. For each 
observation, First, we extract the 4 ms time bin res- 
olution light curve of Channel 0-35 (corresponding to 
2-11.2 keV, 2-13 keV for Epoch II and Epoch III re- 
spectively). Second, we divide the light curve into 32s 
segments and sort them by their count rate. Third, 
we average these segments into several flux bins in or- 
der of the count rate with each flux bin contains ~ 20 
segments (Actually the rms-flux relation does not vary 
with the time-scale of the flux bin, this is also consis- 
tent with the result in Cyg X-l that the linear rms- 
flux relation is valid not only on short time-scales b ut 
also on longer time-scales, see Gleissner et al. ( 20041 )). 
Fourth, we produce the PDSs of flux bins using the 
program 'POWSPEC 1.0' with the Poisson noise sub- 
tracted and the power normaliz ed to squared mean in- 
tensity ( ((Miyamoto et al.lll991l) ). Finally, we integrate 
the power in a certain frequency range and multiply 
the square root of the integral value by the flux (count 
rate) of the flux bin, i.e., 



1/2 



(P(f))df 



F. 



(1) 



where a means the absolute rms, (P(f)) is the av- 
erage power of a flux bin, which can read out from the 
PDS, and F is the count rate of the flux bin (normal- 
ized to one PCU), f± and ]i are the lower and upper 
integral limits fixed at 0.125 Hz and 60 Hz respectively 
so as to include the contribution of QPOs. In this way, 
we can get the rms and flux of all the flux bins of an 
observation and thus derive the rms-flux relation. 

The error of rms is given by the following equation: 



(Aa) 2 = ( 



da 



da. 



WU-/ {AP{f))2 + i dF )2{AF)2 



(2) 



where AP(/) is error of P(f), and AF is the error of the 
flux. As us ual, we use the count rate of each flux bin 
as th e flux ( duttlev and McHardvlboOll lGleissner et al 



2004 iLi et ailboioHHeil and VaugharJboioh ). 

In this paper we have not corrected the effect of dead 
time of the instruments, because it accounts for only 
~1.2% of the total rms and do esn't affect our resu lts, 
(we use the dead time model of Z hang et al. (|l995l ) to 



compute the instru mental dead t i me m odification. See 
also equation (4) in iNowak et al.l (|1999l ).) 



3.1 rms-flux relations in Proposal P10066 

As a example, we will show the rms-flux relation of Pro- 
posal P10066 in this section, because it clearly shows 
three branches in the CCD. Proposal P10066 includes 
three observations, i.e. 10066-01-01-000, 10066-01-01- 
001 and 10066-01-01-00. We display its CCD and hard 
color-intensity diagram (HID) in Fig. [TJ The CCD 
shows three branches clearly, i.e., the top, middle and 
bottom branch. Their PDS are shown in Fig. [2] Be- 
cause both the top and middle branches exhibit QPOs 
between ~15 to ^60 Hz, they locate in the HB. Given 
the shape of the top branch, we name it VHB, and call 
the middle branch HB. The bottom branch is NB bc- 
cause of a power law PD S and its ^7 Hz QPO (see 
Hasi nger and van der Klisl ril989h ). Using the method 



in Section 2.2, we obtain the rms-flux relations of the 
three branches, and the results are shown in Fig. [3l 

Our analysis shows that the rms is negatively cor- 
related with flux in the VHB, but positively correlated 
in the NB of P10066 (Fig. [3). In the HB of P10066, 
the rms-flux relation is not a linear correlation any- 
more. In fact, the rms positively correlated with flux 
in the left part of HB and negatively correlated in the 
right. In this manner the rms has a maximum value 
at the flux of ^1400 cts/s, approximately correspond- 
ing to the middle part of HB. For comparison we show 
the whole rms-flux relation of Proposal P10066 in Fig. 
[U By using a linear model to fit the relation of each 
branch in Fig. |4l we derive the turn-over points be- 
tween different branches. In the VHB of PI 0066, the 
flux is less than ^950 cts/s/pcu and the QPO frequency 
of each flux bin is less than ^30 Hz. In the left HB of 
P10066, the flux is more than ^950 cts/pcu and less 
than ~1430 cts/s/pcu, and its QPO frequency is be- 
tween -30 Hz and -50 Hz. The flux of the right HB 
of P10066 is more than — 1430cts/s/pcu and its QPO 
frequencies range from ~50 Hz to ~55 Hz. Those turn- 
over points are listed in Table [2] 

3.2 rms-flux relations in the other observations 

We derive the rms-flux relations of all the other obser- 
vations listed in Table [TJ The results are shown in Fig. 
EED Fig. show the results of ObsID 30046-01-12- 
00, who locates in the VHB of the CCD. Its flux bin 
frequencies and flux are respectively less than ^30 Hz 
and ~950cts/s/pcu, and the rms-flux relation is neg- 
ative. These results are similar to those of the VHB 
P10066. Fig. |5] shows the results of these observations 
locating in HB. There rms-flux relations are the same as 
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in the HB of P10066. By fitting their rms-flux relations, 
we derive the turn-over points are ^1460cts/s/pcu, 
~1480cts/s/pcu, ~1430cts/s/pcu and — 1400cts/s/pcu 
for 10065-02-01-000, 10065-02-01-001, 10065-02-01-002, 
and 10065-02-01-003, respectively. Moreover, their 
QPO frequencies are between ^38 Hz and ^50 Hz for 
the left HB, and between ~50 Hz and ~55 Hz for the 
right HB (Table EJ). The ObsIDs 20053-04-01-010 and 
20053-04-01-020 locate in the NB and their results are 
shown in Fig [71 There two observations do not have 
obvious QPOs but a bump at ^5-7 Hz, and their hard 
colors are less than ~0.45. Although their flux, between 
^950cts/s/pcu and ^1100cts/s/pcu, are less than in 
the NB of P10066, their rms-flux relations are positive. 

Fig. [8] shows the results of these observations locate 
only in the left part of HB. For the ObsID 30046-01- 
03-00, its flux is less than ~1300cts/s/pcu and its fre- 
quency is between ~45 Hz and ~53 Hz. However, the 
situation is a little different for ObsID 300418-01-01- 
00. Its flux is less than ^1500cts/s/pcu and the QPO 
frequencies range from ^30 Hz to ~48 Hz. As a re- 
sult, their rms-flux relations are positive. There are 
three observations locate in the right part of HB(see 
Table Q]) and their results are shown in Fig. [9] Their 
hard colors are more than ^0.45 while their flux are be- 
tween ~800cts/s/pcuand ^1100cts/s/pcu. For ObsIDs 
20053-04-01-03 and 20053-04-01-04, the QPO frequen- 
cies are around ^53-55 Hz, but it locates around ~45 
Hz in the ObsID 20053-04-01-030. The same as in the 
right HB of P10066, their rms are negatively correlated 
with flux. Because these observations shown in Fig. [8} 
[9] occupy only a portion of HB, their rms-flux relations 
show only a portion of the arch-like relation of the HB 
of P10066 either. 

In addition, there is one observation, 30046-01-08- 
00, locating at the vertex of the VHB and HB (see Fig. 
[T0|) . This observation has two branches. The flux of the 
VHB part is less than ^960 cts/s/pcu, and its QPO 
frequencies range from ^26 Hz to ^35 Hz, while the 
flux of the left HB part is more than ~ 960 cts/s/pcu 
and its QPO frequencies are between ^37 Hz and ^50 
Hz(see Tabled]). As shown in the right panel of Fig. 
[lOj the VHB part shows a negative correlation in the 
rms-flux plot, and the left HB part shows a positive 
one, which are consistent with those results above. 

3.3 the correlations between rms and spectral 
parameters in P10066 

Because the other observations have the same results 
of rms-flux relations as in Proposal P10066, our spec- 
tral fitting concentrate only on Proposal PI 0066. First, 
we split the data into 25 intervals according to its po- 
sition on the CCD. The average hard color (HC) and 



soft color (SC) of the 25 intervals are listed in Table 
[3l (Because there is no need to correct the gain dif- 
ference between Epoches II and III, we haven't cor- 
rected the colors in this subsection and they are differ- 
ent from above.) Second, we compute their rms using 
the method described above (see Table [3]). Finaly, we 
extract their spectra with the Standard 2 mode data 
and use the X-ray spectral-fitting program, XSPEC, to 
analyze their spectra. 

Considering the controversy over the location and 
nature of the X-ray emission regions, we fit their 
spectra with the Eastern and Western model, re- 
spectively. We use the absorbed disk blackbody, 
Comptonization plus an extra Gaussian line (wabs( 
diskbb+compTT-|-gauss)) as the Eastern model, and 
the absorbed blackbody, Comptonization plus an ex- 
tra Gaussian line (wabs(bb+compTT+ga uss)) as the 



Weste rn model, i.e. the Birmingham model lChurch and Balucinska-Chi 
I l995h . The fitting energy band is 3-30keV and the ab- 
sorption column density N^ is fixed at 0.2 xlO 22 . Ac- 
cording to their reduced x 2 values shown in Table El we 
can find that the Western model gives a better fitting 
results for those intervals whose ID > 17. We also give 
the inner radius of the accretion disk of the Eastern 
Model in Table |3] by a ssuming an inclinati on angle of 
the system of 60° (see (|Di Salvo et al.ll2002i ) V Because 
the inner disk radius is less than the typical radius of 
a NS, i.e., ^10km, the Eastern Model isn't reasonable 
for Cyg X-2 in this situation. The spectral parameters 
of the Birmingham model are listed in Table HI The 
relations between the rms and flux of these two com- 
ponents of Birmingham model is shown in Fig. 111! We 
also present the relation between the rms and flux ratio 
of blackbody and Comptonization components in Fig. 

m 



4 Discussion 



Li et al. I (|2010h firstly reported a negative rms-flux re- 
lation in Cyg X-2, and then it is als o found in the bla ck 
hole XRB XTE J1550-564 (see also lHeil et ail |201lh h 
The positive relation can, in principle, be readily un- 
derstood: the larger the flux, the higher the accretion 
rate which leads to a higher level of instability. Never- 
theless, the negative correlation between rms and flux 
in Cyg X-2 a nd XTE J 1550- 56 4 is hard t o fit in to this 
scenario (see iLi et all (|2010l ); iHeil et all (|201ll )). As 
shown in Fig. [TU1 the complicated rms-flux relation of 
Cyg X-2 leads to a maximal rms in the middle of HB, 
which means that Cyg X-2 has a maximum level of in- 
stability no matter how luminous it is. Moreover, the 
negative relation dose not change with the integral fre- 
quency range. In other words, if we do not include the 
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QPO and only integral from 0.125 Hz to 10 Hz when we 
compute the rms, the rms-flux relation is also negative 
in the VHB. This means that the negative relation is 
not caused by the QPO. We also normalize the rms by 
dividing the flux of each flux bin, and the normalized 
rms shows a positive correlation with flux in the NB, 
but negative in the VHB, which is totally the same as 
not normalized. The same results mean that the rms- 
flux relation isn't due to the flux change of different flux 
bins. 

To compare NSXRB Cyg X-2 with BHXRB Cyg 
X-l, we investigate the rms-flux relations in differ- 
ent branches in Cyg X-2. As shown above, the rms 
is positively correlated with flux when Cyg X-2 lo- 
cates in the NB or at the left part of HB. However, 
it is negatively correlated with flux in the VHB or at 



the right 
ported in 



part of HB. This is different from those re- 



Uttlcv and McHardyj (120011); lUttlevI (120 04); 



Heil and Vaughan ( 2010h . Especially, Gleissner et al 



( 20041 ) found that the positive rms-flux relation is valid 
throughout all states of the BHXRB Cyg X-l. More- 
over, the slope of the rms-flux relation of the hard state 
is steeper than in the soft and intermediate states of 
Cyg X-l. However, the slopes of the hard state (VHB 
state and its slope is negative) is flatter than in the soft 
state (NB state, its slope is positive) in Cyg X-2. We 
guess that this difference may come from the compact 
star. 

However, although Cyg X-l and Cyg X-2 have dif- 
ferent results in the hard and soft states, the rms-flux 
relation shows some similarities if we assume HB is the 



intermediate state of Cyg X-2. Ideissner et al.l (|2004l) 
reported an 'arch'-like rms-flux relation in the interme- 
diate state of Cyg X-l, but they claimed that the 'arch' 
relation is due to the buffer overflow. 

The so-called buffer overflow occurs when a program 
or process tries to store more data in a buffer (tem- 
porary data storage area) than it could hold. Due to 
the softness of the X-ray spectrum, the likelihood for 
such overflows is the highest in the lowest energy band. 
Still, for the different configuration, the possibility of 
overflow is lower when the configuration threshold is 
higher. In this work, we also find some 'arch' relations 
in the intermediate state of Cyg X-2 and try to deter- 
mine wether the arch shape in Cyg X-2 is due to the 
buffer overflow or not. We've tried the Event Encoded 
mode data (see Table Q] and Fig. [6]), the higher thresh- 
old mode data, and higher energy band data, which 
are less likely suffered from buffer overflows (also see 



http://heasarc.gsfc.nasa.gov/docs/xte/RXTE_tech_apl 



pend.pdf)), and found that these data also show simi- 
lar arch relations. On the other hand, we've done some 
simulations and the simulation results don't support 



the argue that the arch relation is due to the buffer 
overflow. As a result, we conclude that the 'arch' rms- 
flux relation may be the characteristic of the state tran- 
sition, which always take place at the QPO frequency 
~50 Hz (see Table 0). 

In the 1980s, two models were proposed to ex- 
plain the X-ray emission of the NSXRB: the Western 
model and the Eastern model. The Eastern model con- 
tains a multi-temperature blackbody from the inner 
disc, plus a Comptonized emissio n of the seed photon s 



from the blackbody of the NS (|Mitsuda et al.l Il989l) 



while the Western model contains a Comptonization 
component and sometimes a second one (i.e., black- 



body emission) for some luminous sources (jWhite et al 
19881) . In the 1990s, a new model - the Birming- 



ham model, closely related to the Western model - is 
proposed. This model assumes that all LMXBs have 
two continuum components: blackbody emission from 
the NS surface and Comptonized emission from an ex- 
tended accretion disc corona above th e accretion disc 
(jChurch and Balucinska-Churchl Il995l ) . Our spectral 
results show that the Birmingham model is more rea- 
sonable than the Eastern model (see Table _J|) , and this 
suggests that the two components in Cyg X-2 are most 
likely the blackbody coming from the surface of the NS 
and the Comptonization generated from the corona. 

From Table _H we know that the temperatures of the 
seed photons, the blackbody and the corona, the optical 
depth, and the flux of the Comptonization component 
change little, while the blackbody flux increases appar- 
ently in the VHB. A scenario hinted from these results 
may be that the corona is not efficiently cooled by the 
soft disk emission, and the latter is also likely stable 
because of the little change in TO. As a result most ma- 
terial move onto the surface of NS, and the rms changes 
obviously. In other words, the negative rms-flux rela- 
tion of the VHB is supposed to contribute from the 
blackbody component (Fig. [TTj) . Similarly, the rms is 
dominated by the Comptonization component in the 
NB (Fig. [TT]). To sum up, the blackbody component 
dominates the rms variability in VHB and its contribu- 
tion to rms decreases as Cyg X-2 evolving through the 
HB to NB, whereas the rms contributed by the Comp- 
tonization component increases and dominates the rms 
variability in NB at last. This inference is also sup- 
ported by the results that the rms negatively correlated 
with flux in energy bands 2-6keV a nd 6-10keV, an d 
positively correlated in 10-21keV (see Li et al. ( 20101) ). 
Because of these totally different manners of these two 
components, Cyg X-2 has a complex rms-flux relation 
along with the CCD track. 

Fig. [___ shows that the instability of Cyg X-2 is 
related to the flux ratio between the blackbody and 
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Comptonization components. Only when those two 
components reach a certain ratio, Cyg X-2 reaches its 
maximum instability. This can be explained in the 
Birmingham model. In this model, the radiation of 
the blackbody shines on the accretion disc and forms 
a corona. On the one hand, the rms produced by the 
blackbody decreases when the flux of the blackbody 
increases. On the other hand, the radiation of the in- 
creasing blackbody will expand the corona region and 
increase the flux of the Comptonization component, 
which will finally increase its rms. As a result, the rms 
has a maximum only when their flux reaches a certain 
ratio. 
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The CCD (left panel) and HID (right panel) of 
Proposal P 10066. The hard and soft color is defined as the 
Nowak, M.A., Vaughan, B.A., Wilms, J., Dove, J.B., Begel- flux ratio between 11.2-18.8keV and 7.5-11. 2kev, and be- 
man, M.C.: Astrophys. J. 510, 874 (1999). doi:10.1086/306610fveen 4.7-7.5keV and 2-4.7keV, respectively. The Intensity 

is defined as the count rate of 2-18. 8keV. 
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Table 1 The information of those observations used in this work. 



ObsID 


date 


epoch duration 


data mode 




Branch 


10065-02-01-000 


23/3/96 


II 


14ks 


E.125us_64M_0.1s 


HB 


10065-02-01-001 


23/3/96 


II 


25ks 


E.125us_64M_0.1s 


HB 


10065-02-01-002 


23/3/96 


II 


19ks 


B.2ms_16A_0_35. 


-Q 


HB 


10065-02-01-003 


23/3/96 


II 


15ks 


B.2ms_16A_0_35. 


-Q 


HB 


10066-01-01-00 


27/3/96 


II 


16ks 


B.2ms_16A_0_35. 


-Q 


VHB 


10066-01-01-000 


26/3/96 


II 


16ks 


B.2ms_16A_0_35. 


-Q 


NB 


10066-01-01-001 


27/3/96 


II 


18ks 


B.2ms_16A_0_35. 


-Q 


HB 


20053-04-01-010 


1/7/97 


III 


25ks 


B_2ms_16A_0_35. 


-Q 


NB 


20053-04-01-020 


1/7/97 


III 


15ks 


B.2ms_16A_0_35. 


-Q 


NB 


20053-04-01-03 


2/7/97 


III 


12ks 


B_2ms_16A_0_35. 


-Q 


right HB 


20053-04-01-030 


2/7/97 


III 


13ks 


B_2ms_16A_0_35. 


-Q 


right HB 


20053-04-01-04 


2/7/97 


III 


15ks 


B.2ms_16A_0_35. 


-Q 


right HB 


30046-01-03-00 


28/7/98 


III 


15ks 


B.2ms_16A_0_35. 


-Q 


left HB 


30046-01-08-00 


30/8/98 


III 


12ks 


B_2ms_16A_0_35. 


_Q VHB-HB vertex 


30046-01-12-00 


25/9/98 


III 


20ks 


B.2ms_16A_0_35. 


-Q 


VHB 


30418-01-01-00 


2/7/98 


III 


12ks 


B_2ms_16A_0_35. 


-Q 


left HB 



Table 2 The turn-over points between different branches. R and f stand for the flux in the unit of cts/s/pcu and the 
frequency of the QPO or the bump in the unit of Hz, respectively. 



ObsID 



VHB 



ft HB 



ht HB 



NB 



10065-02-01-000 
10065-02-01-001 
10065-02-01-002 
10065-02-01-003 
VHB of P10066 
left HB of PI 0066 
right HB of PI 0066 
NB of P10066 
20053-04-01-010 
20053-04-01-020 
20053-04-01-03 
20053-04-01-030 
20053-04-01-04 
30046-01-03-00 
30046-01-08-00 
30046-01-12-00 
30418-01-01-00 



R<950, f<30 



R<1460, 38<f<50 
R<1480, 38<f<50 
R<1430, 38<f<50 
R<1400, 38<f<50 

950<R<1430, 30<f<50 



R<960, 26<f<35 
R<950, f<30 



R<1300, 45<f<53 
960<R, 37<f<50 

R<1500, 30<f<48 



1460<R, 50<f<55 
1480<R, 50<f<55 
1430<R, 50<f<55 
1400<R, 50<f<55 



1430<R, 50<f<55 



5<f<7 
5<f<7 
5<f<7 



800<R<1100, f~53-55 
800<R<1100, f~53-55 
800<R<1100, f~45 
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Table 4 The fitting parameters of the Birmingham model. The fitting energy band is 3-30keV and the absorption column 
density N h is fixed at 0.2 xlO 22 . 



ID 




TO 1 






Tbb 1 






kT 1 






T 2 






Fe 3 




flux& 


^ UX iompTT 


fluxf ^ , 

total 


1 





183±0 


005 


1 


39±0 


07 


3 


11±0 


05 


5 


67±0 


16 


6 


44±0 


09 


i 


uu 


R 
u 




7.84±0.01 


2 





194±0 


004 


1 


27±0 


05 


3 


02±0 


04 


5 


94±0 


12 


6 


46±0 


09 


n 

u 


QQ 

■j -j 


u 


81 


7.96±0.01 


3 





202±0 


004 


1 


26±0 


05 


3 


03±0 


04 


5 


91±0 


12 


6 


49±0 


08 


I 


07 


u 


77 


8.00±0.01 


4 





186±0 


004 


1 


26±0 


05 


3 


02±0 


04 


5 


89±0 


13 


6 


58±0 


08 


I 


1 3 




76 


8.02±0.01 


5 





254±0 


003 


1 


26±0 


04 


3 


01±0 


04 


5 


90±0 


12 


6 


50±0 


08 


I 


1 3 




80 


8.07±0.01 


6 





202±0 


004 


1 


29±0 


05 


3 


02±0 


05 


5 


80±0 


14 


6 


48±0 


08 


I 


1 7 


P. 

\J 


84 


8.13±0.01 


7 





253±0 


003 


1 


21±0 


03 


2 


96±0 


03 


5 


91±0 


11 


6 


56±0 


09 


I 


32 


u 


78 


8.23±0.01 


8 





215±0 


003 


1 


19±0 


03 


2 


94±0 


03 


5 


89±0 


10 


6 


49±0 


09 


I 


33 


u 


82 


8.28±0.01 


9 





178±0 


002 


1 


17±0 


03 


2 


93±0 


02 


5 


95±0 


10 


6 


59±0 


08 


1 

_L 


to 


a 
\j 


ou 


8.35±0.01 


10 





214±0 


003 


1 


21±0 


03 


2 


95±0 


03 


5 


77±0 


11 


6 


50±0 


10 


I 


BO 


a 
u 


ou 


8.57±0.01 


11 





192±0 


004 


1 


18±0 


02 


2 


85±0 


04 


G 


02±0 


14 


6 


56±0 


11 


o 


oo 


7 
i 




9.90±0.01 


12 





177±0 


004 


1 


17±0 


02 


2 


85±0 


04 


G 


04±0 


13 


6 


59±0 


14 


2 


75 


7 


42 


10.27±0.01 


13 





184±0 


003 


1 


11±0 


02 


2 


73±0 


03 


G 


47±0 


12 


6 


53±0 


16 


2 


71 


7 


53 


10.37±0.02 


14 





219±0 


003 


1 


14±0 


02 


2 


78±0 


03 


G 


31±0 


12 


6 


62±0 


12 


2 


85 


7 


57 


10.53±0.02 


15 





228±0 


003 


1 


14±0 


02 


2 


76±0 


03 


G 


37±0 


12 


6 


57±0 


14 


2 


88 


7 


77 


10.77±0.01 


16 





185±0 


003 


1 


12±0 


01 


2 


73±0 


03 


G 


60±0 


12 


6 


60±0 


14 


3 


06 


7 


81 


11.01±0.01 


17 





204±0 


004 


1 


13±0 


02 


2 


70±0 


03 


G 


72±0 


13 


6 


33±0 


18 


3 


40 


8 


41 


11.98±0.02 


18 





207±0 


005 


1 


17±0 


02 


2 


70±0 


03 


G 


76±0 


15 


6 


58±0 


15 


4 


23 


8 


91 


13.22±0.02 


19 





182±0 


005 


1 


21±0 


02 


2 


75±0 


04 


G 


44±0 


16 


6 


64±0 


16 


4 


80 


9 


28 


14.16±0.02 


20 





184±0 


005 


1 


18±0 


01 


2 


68±0 


03 


G 


86±0 


15 


6 


59±0 


17 


5 


03 


9 


59 


14.70±0.02 


21 





181±0 


004 


1 


19±0 


01 


2 


69±0 


03 


G 


33±0 


11 


6 


64±0 


18 


5 


76 


8 


97 


14.80±0.02 


22 





190±0 


004 


1 


16±0 


01 


2 


63±0 


03 


G 


24±0 


11 


6 


62±0 


11 


5 


42 


8 


23 


13.74±0.02 


23 





184±0 


003 


1 


12±0 


01 


2 


56±0 


03 


G 


02±0 


11 


6 


62±0 


13 


5 


28 


7 


07 


12.46±0.02 


24 





196±0 


003 


1 


07±0 


01 


2 


41±0 


02 


G 


03±0 


10 


6 


57±0 


09 


5 


07 


6 


12 


11.32±0.01 


25 





214±0 


004 


1 


00±0 


01 


2 


26±0 


02 


G 


49±0 


10 


6 


30±0 


07 


4 


69 


5 


59 


10.55±0.02 



Their units are keV and represent the temperature of seed photons, blackbody and the corona respectively. 

2 r stands for the optical depth. 

3 Fe denotes the energy of the iron line. 

4 Their units are 1CF 9 ergss~ 1 cmT 2 and represent the 3-30keV flux of blackbody, Comptonization component and 
total flux respectively. 




0.01 0.1 



Fig. 2 The PDSs of VHB(left), HB(middle) and NB(right) for Proposal P10066. 
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Fig. 3 The rms-flux relations of VHB(left), HB(middle) and NB(right) for Proposal P10066 
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Fig. 5 This figure shows the CCD(left), PDS(middle) and rms-flux relation (right) of ObsID 30046-01-12-00, which locates 
at the VHB and displays a negative correlation. 
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Fig. 4 The complete rms-flux relation of Proposal P10066. 
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Fig. 11 The rms-flux plot. The X-axis of the left panel 
is the flux of blackbody, and the right is the flux of Comp- 
tonization component. 
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Fig. 12 The rms-ratio plot of Proposal P10066. The X- 
axis stands for its flux ratio between blackbody and Comp- 
tonization components. 
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Fig. 6 The CCDs (left column), PDSs (middle column) and rms-flux relations (right column) of those observations 
locating at the HB. Row 1, 2, 3 and 4, show the relults of ObsIDs 10065-02-01-000, 10065-02-01-001, 10065-02-01-002 and 
10065-02-01-003, respectively. The left HB is marked with triangles and the right HB with squares in each CCD. 
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Fig. 7 The CCDs (left column), PDSs (middle column) and rms-flux relations (right column) of observations locate in 
the NB. The top and bottom panels show the results of ObsIDs 20053-04-01-010 and 20053-04-01-010 respectively. 
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Fig. 8 The CCDs (left column), PDSs (middle column) and rms-flux relations (right column) of observations locate in the 
left part of HB. The top and bottom panels show the results of ObsIDs 30046-01-03-00 and 30418-01-01-00 respectively. 
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Soft Color Frequency (Hz) Flux (counts s~ 1 pcir 1 ) 

Fig. 9 The CCDs (left column), PDSs (middle column) and rms-flux relations (right column) of observations locate in 
the right part of HB. The top, middle and bottom panels show the results of ObsIDs 20053-04-01-03, 20053-04-01-030 and 
20053-04-01-04, respectively. 
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Fig. 10 The CCD (left), PDS (middle) and rms-flux relation (right) of ObsID 30046-01-08-00, which locates in the vertex 
of VHB and HB. 



14 



Table 3 The \ °f the Western and Eastern models. Ri„ 
is the inner disk radius derived from the Eastern Model by 
assuming an inclination angle of the system of 60°. 



ID 


HC 


SC 


Western 1 


Eastern 2 


rms 




1 


U.oyo 





885 


u.y/(oo j 


U.yl(i3D ) 


yi.o±u.4 


Q Q 
o.o 


n 
Z 


n QO/i 
U.oy4 





885 


U.oo(oo ) 


a on/ 
U.oy(OD ) 


yi.o±u.4 


A O 


o 

d 


A OA1 

0.391 





883 


0.89(56) 


A A 1 i £ O \ 

0.91(56) 


AA A I A A 

90.0±0.4 


4.6 


A 

4 


U.ooy 





875 


l.lo(DO ) 


l.z4(0D ) 


qa £j_n A 


A Q 

4.o 


r 



n qca 





876 


l.Uo(oo ) 


l.U< (00 ) 


qq £_i_n a 


A V 

4. i 


a 
u 


u.ooo 





872 


U.OU^OD ) 




ni K_i_n 7 
yi.oznu. i 


A A 


7 


0.377 





856 


0.55(56) 


0.56(56) 


86.2±0.4 


5.8 


8 


0.373 





850 


0.62(56) 


0.65(56) 


83.4±0.4 


6.1 


9 


0.372 





848 


0.96(56) 


0.99(56) 


85.0±0.4 


6.6 


10 


0.364 





839 


1.00(56) 


1.04(56) 


82.8±0.6 


6.3 


11 


0.343 





830 


0.82(56) 


0.86(56) 


85.9±1.0 


8.3 


12 


0.340 





832 


0.57(56) 


0.65(56) 


87.7±0.7 


8.5 


13 


0.339 





830 


0.74(56) 


0.86(56) 


89.1±0.7 


9.1 


14 


0.339 





833 


0.69(56) 


0.83(56) 


90.7±0.7 


8.9 


15 


0.338 





839 


0.65(56) 


0.77(56) 


91.2±0.7 


9.0 


16 


0.338 





840 


0.76(56) 


1.11(57) 


93.0±0.7 


8.9 


17 


0.334 





850 


0.73(56) 


1.41(57) 


100.9±1.4 


8.8 


18 


0.329 





857 


0.98(56) 


1.55(57) 


108.4±1.1 


8.9 


19 


0.321 





858 


0.62(56) 


1.37(57) 


106.7±0.9 


8.9 


20 


0.322 





861 


0.47(56) 


1.42(57) 


104.7±0.8 


9.1 


21 


0.295 





835 


0.73(56) 


1.57(57) 


69.4±0.6 


10.6 


22 


0.282 





812 


0.62(56) 


1.61(57) 


56.7±0.7 


11.3 


23 


0.258 





768 


0.77(56) 


1.72(57) 


42.0±0.9 


12.4 


24 


0.233 





726 


1.06(56) 


2.27(57) 


25.6±1.4 


13.9 


25 


0.220 





702 


1.61(56) 


3.63(57) 


20.9±1.2 


15.2 



Western Model, i.e., wabs(bb+compTT+gauss). 

2 Eastern Model, i.e., wabs(diskbb+compTT+gauss). 

3 Rin is derived from the Eastern Model in unit of km. 



